Surface polaritons on left-handed cylinders: 
A complex angular momentum analysis 
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We consider the scattering of electromagnetic waves by a left-handed cylinder - i.e., by a cylinder 
fabricated from a left-handed material - in the framework of complex angular momentum techniques. 
We discuss both the TE and TM theories. We emphasize more particularly the resonant aspects 
of the problem linked to the existence of surface polaritons. We prove that the long-lived resonant 
modes can be classified into distinct families, each family being generated by one surface polariton 
propagating close to the cylinder surface and we physically describe all the surface polaritons by 
providing, for each one, its dispersion relation and its damping. This can be realized by noting 
that each surface polariton corresponds to a particular Regge pole of the S matrix of the cylinder. 
Moreover, for both polarizations, we find that there exists a particular surface polariton which 
corresponds, in the large-radius limit, to the surface polariton which is supported by the plane 
interface. There exists also an infinite family of surface polaritons of whispering-gallery type which 
have no analogs in the plane interface case and which are specific to left-handed materials. 

PACS numbers: 78.20.Ci, 41.20.Jb, 73.20.Mf, 42. 25. Fx 



I. INTRODUCTION 

In an article published in 1967 (see Ref. Q for the En- 
glish translation), Veselago speculated upon the alter- 
ation of electromagnetic and optical phenomena in hy- 
pothetical materials for which the electric permittivity 
and the magnetic permeability (and therefore the refrac- 
tive indesi^) were simultaneously negative in a certain 
range of frequencies. He predicted more particularly the 
existence, in such materials, of anomalous effects such 
as a reversed Doppler shift, reversed Cerenkov radiation, 
negative radiation pressure, and inverse Snell-Descartes 
law. Since in such substances the electric field, the mag- 
netic field and the wave vector of a plane wave form a 
left-handed system, so that the Poynting vector and the 
wave vector have opposite direction, he referred to them 
as left-handed media and here we shall use this termi- 
nology, even if today other authors prefer to call this 
kind of materials "negative index media" , "double nega- 
tive media" or "Veselago's media". Four years ago, fol- 
lowing insights from Pendry and coworkers^^, Schultz, 
Smith and colleaguesSi^ have been able to build, for the 
first time, an artificial left-handed medium by combin- 
ing arrays of wires and split-ring resonators and to ex- 
perimentally test "left-handed electromagnetism" in the 
microwave frequency range. Since then, several other 
groups have successfully fabricated left-handed media 
and experimentally studied the alteration of electromag- 
netic phenomena (see, for a comprehensive list of refer- 
ences, the popular article by Pendry and SmithS). Now, 
the possibility to create left-handed materials working 
in the optical domain is seriously foreseen and, very re- 
cently, negative refraction has been even observed at 
infrared wavelengths in the context of photonic crys- 



tal physicsiS. Of course, the unusual and remarkable 
properties of left-handed media could revolutionize op- 
toelectronics and communications. Many technological 
applications are already considered including superlenses, 
bandpass filters, beam guiders and light-emitting devices. 
As a consequence, this recent new field of physics has at- 
tracted the interest of many researchers and is rapidly 
evolving and the corresponding scientific literature is ex- 
ploding. 

In this article, we shall focus our attention on a particu- 
lar problem of left-handed electromagnetism: namely, the 
scattering of an electromagnetic wave by a "left-handed 
cylinder" - i.e., by an homogeneous cylinder fabricated 
from a left-handed material. This problem has been al- 
ready considered by Kuzmiak and Maradudinii but we 
shall reexamine it in the framework of complex angular 
momentum (CAM) techniquesiSiii In fact, we are above 
all interested in the resonant aspects of scattering linked 
to the existence of surface polaritons (SP's) and, by using 
CAM techniques, we shall be able to provide a theoretical 
description of SP's propagating close to the left-handed 
cylinder surface as well as a physical explanation for the 
excitation of the associated resonant modes. 

Let us recall here that SP's are electromagnetic sur- 
face waves propagating close to the interface separating 
two different media with an amplitude that decays in 
an exponential fashion perpendicularly to the interface 
and into both media. SP's supported by metal-dielectric 
or semiconducting-dielectric interfaces have been exten- 
sively studied during the last 40 years because of the fun- 
damental role they play in the context of the interaction 
of electromagnetic radiation with matter but also because 
of their numerous practical applications in physics, chem- 
istry and biology (for reviews on this subject we refer to 
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Refs. HUll Il6lll7fl . Recently, activity has focused par- 
ticularly on the role of SP's in pho to nic crystal physics 
(see, for example, Refs. [TMl9ll2nll2lll^l23l24l25l2fl27l2^ 

as well as on their role in the enhanced transmission of 
light through periodic arrays of holes in a metal film 29 
(see also Ref. |21J and references therein). Of course, the 
scientific literature dealing with SP's localized at the in- 
terface separating a left-handed medium and a conven- 
tional one is much less developed and it is still too early to 
judge their importance in left-handed electromagnetism. 
However, it is already obvious that it will not be negli- 
gible: indeed, it seems that SP's play a central role in 
the superlensing phenomenon2ii22i2&i2i and in the giant 
Goos-Hanchen effect recently discovered 35 . 

The theory of SP's suppor ted by le ft-handed me- 
dia has been studied in Refs. I36l37l38l in the case of 
a plane interface. In that configuration, we can con- 
sider that the properties of SP's are now completely 
known. They have been obtained from rather elementary 
calculations involving homogeneous and inhomogeneous 
plane waves. The existence of SP's in the presence of a 
curved interface has been noted in the article by Kuzmiak 
and Maradudinil previously cited and in an article by 
RuppmSS, dealing with scattering from a sphere. But it 
should be noted that neither of these two articles provide 
a clear physical description of SP's. This is mainly due 
to the fact that, unlike the case for a plane interface, in 
the case of a curved interface the wave equation cannot 
be solved in terms of elementary functions. As a conse- 
quence, the description of SP's involves transcendental 
equations which obscures any analysis. Fortunately, it is 
possible to circumvent these difficulties by using CAM 
techniques (here, we refer to the Introduction of Ref. 
for a short bibliography on this topic) in connection with 
modern aspects of asymptotics^i^2i2ii2i. It is exactly 
what we did in a recent article^ dealing with scattering 
of electromagnetic waves from metallic and semiconduct- 
ing cylinders - i.e., in the presence of dispersive materials. 
This allowed us to provide a clear physical explanation 
for the excitation mechanism of the resonant surface po- 
lariton modes (RSPM's) as well as a simple mathematical 
description of the SP's that generate them. In the present 
paper, we shall extend this analysis to left-handed cylin- 
ders but before we embark on this work we shall briefly 
discuss the limitations of our approach. We first recall 
that the CAM method is an asymptotic approach and 
that formally it is only valid for high frequencies or more 
precisely when the radius a of the cylinder is large com- 
pared to the wavelength A = 2ttc/oj of the electromag- 
netic field. In fact, this limitation must not be taken too 
seriously. Indeed, in practical applications CAM tech- 
niques provide good results even if A ~ o. Moreover, 
in the present problem we cannot assume that a is very 
large compared to A. Indeed, in that case, the wavelength 
of the electromagnetic field could become comparable to 
the size a' of the unit cell of the left-handed material 
and the cylinder could not be treated as homogeneous. 
In other words, the analysis we shall provide in the next 



sections is formally valid for a' <C A < a. 

Our paper is organized as follows. In Sec. II, we intro- 
duce our notations and we construct the S matrix of the 
system. We consider the TE theory (H polarization) as 
well as the TM theory (E polarization). We then discuss 
the resonant aspects of our problem for both theories. 
In Sec. Ill, by working in the CAM plane, we quali- 
tatively describe the SP's supported by the left-handed 
cylinder. Then, in Sec. IV, by using CAM techniques, 
we establish the connection between these SP's and the 
associated RSPM's. In Sec. V, by using asymptotic tech- 
niques, we describe semiclassically the different SP's and 
we provide analytic expressions for their dispersion re- 
lations and their damping. We show more particularly 
the existence of SP's of whispering-gallery type. Finally, 
in Sec. VI, we conclude our article by emphasizing the 
main results of our work. 



II. EXACT S MATRICES AND SCATTERING 
RESONANCES FOR THE TE AND THE TM 
THEORIES 

From now on, we consider the scattering of electromag- 
netic waves by a cylinder with circular cross section and 
radius a having an effective frequency-dependent permit- 
tivity e(w) and an effective frequency-dependent perme- 
ability p{uj). Here and in the following, we implicitly 
assume the time dependence exp(— itot) for electric and 
magnetic fields. We consider that the cylinder is embed- 
ded in a host medium of infinite extent having the elec- 
tromagnetic properties of the vacuum. We introduce the 
usual cylindrical polar coordinate system (p,8,z). It is 
chosen so that the cylinder and surrounding medium, re- 
spectively, occupy the regions corresponding to the range 
< p < a (region II) and to the range p > a (region I). 
Furthermore, in order to describe wave propagation we 
also introduce the wave number 

k(u) = - c (1) 

where c denotes the velocity of light in vacuum, and the 
refractive index of the cylinder 

n{ijj) = \J e(w)/x(co>). (2) 

As far as the electric permittivity e(w) and the mag- 
netic permeability /i(w) of the cylinder are concerned, we 
assume they are, respectively, given by 

C ( W ) - 1 - ^ ( 3 ) 

and 

where < F < 1 and Uf, = uj /y/l — F. Here we have 
considered that the cylinder is fabricated from a meta- 
material. Of course, the parameters u> p , loq and F depend 
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on its structure. But we would like to be as general as 
possible in our analysis and therefore we do not restrict 
our study to any particular metamaterial (see, however, 
the discussion in the last paragraph of our conclusion). 
As a consequence, we do not attribute any "microscopic" 
interpretation to the parameters u) p , ujq and F in terms 
of the internal structure of the metamaterial considered. 
Similarly, we do not precise the frequency range where 
u) p and ujq lie. We only assume that loq < u>b < u p . We 
then have e(w) < in the frequency range uj £ ]0,u> p [ 
and n(ui) < in the frequency range u> 6 Jo^t^. Thus, 
the electric permittivity, the magnetic permeability and 
the refractive index are simultaneously negative in the 
region luq < uj < u>b- In that region the metamaterial 
presents left-handed behavior. As far as the numerical 
aspects of our work are concerned, we shall work with 
F = 0.4 and with the reduced frequencies Lo^a/c = 5.52, 
LOba/c rs 7.127 and u p a/c = 11.04. Even though we re- 
strict ourselves to that particular configuration, the re- 
sults we shall obtain numerically are in fact very general 
and they permit us to correctly illustrate the theory. 

Here and from now on, we choose to treat our problem 
in a two-dimensional setting, ignoring the z coordinate. 
We briefly recall the equations governing the TE theory 
where the magnetic field H is parallel to the cylinder axis 
(H polarization) and the TM theory where the electric 
field E is parallel to the cylinder axis (E polarization). 
From Maxwell's equations it is easy to show that the z 
components of the magnetic and electric fields satisfy the 
Hclmholtz equation 



A x + n\uj) (- 
V c 



H?(x) = 

= o 



i?*(x) = 
4(x) = 



for < p < a, 
(5a) 

for p > a, (5b) 



with x = (p, 9) and where the Laplacian A x is expressed 
in the polar coordinate system. For the TE theory, from 
the continuity of the tangential components of the electric 
and magnetic fields - i.e., of Eg and H z - at the interface 
between regions I and II, it can be shown that the z 
component of the magnetic field satisfies, for < 9 < 2tt, 



Hl(p = a,9) = H?(p = a,9), 



dHl 

dp 



(p = a,9) 



i dm 



e{ui) dp 



(ft = a, 



(6a) 
(6b) 



For the TM theory, the continuity of the tangential com- 
ponents of the electric and magnetic fields - i.e., of E z 
and Hg — at the interface between regions I and II per- 
mits us to show that the z component of the electric field 
satisfies, for < 9 < 2tt, 



E I z (p = a,9)=E I z 1 (p = a,6), 



(7a) 



dEl 
dp 



(p = a,9) 



1 dEl 1 



p(u>) dp 



(p = a,9). (7b) 



We can now construct the S matrix for the cylinder for 
both polarizations. Because of the cylindrical symmetry 
of the scatterer, the S matrix is diagonal and its elements 
Sui{w) are given by Sw{uj) = Si(ui) 6w. It should be 
recalled that the 5" matrix is of fundamental importance 
because it contains all the information about the scatter- 
ing process. Its components appear in the Green func- 
tions of the problem, in the scattered field when a plane 
wave excites the cylinder as well as in both the scattering 
amplitude and the total scattering cross section. For the 
TE and the TM theories, we shall denote, respectively, 
by Sf(ui) and Sf(uj) the S matrix diagonal elements. 
For a given angular momentum index t £ 7,, the coef- 
ficients Sf and Sf are respectively obtained from the 
partial wave (H z ) ( and (E z ) e solutions of the following 
problem^: 

(i) (H z ) l and {E z ) l satisfy the Hclmholtz equation 10, 

(ii) (H z )e an d (Ez)ii respectively, satisfy the boundary 

conditions © and Q), 

(iii) at large distance, (H z ) e and {E z ) t respectively 
present the asymptotic behaviors 



(H z ) t (p,0) 



(Ez) t (p,e) 



1 



-i(kp-ln /2-7T/4) 



p^+oo ^/2~nkp 
1 



,-i(fcp-£-7r/2-7r/4) 



+ Sf ( w ) e 'MV2-/4)) e <« 

Outside the cylinder (region I), the solution of Eq. © is 
expressible in terms of Bessel functions (see Ref. l4q) as a 
linear combination of Ji{ujp/ c)e lle and H^\ujp/c)e ie9 . 
Inside the cylinder (region II), it is proportional to 
Ji (n(uj) uip/c)e tle . As a consequence, the partial waves 
{H z ) e and (E z ) e solutions of (i) and (ii) can be obtained 
exactly. Then, by using the standard asymptotic behav- 
ior of Hankel functions (x) and (x) for x — > oo 
(see Ref. l46j) . we find from (iii) the expressions of the di- 
agonal elements Sf and Sf of the S matrix for the TE 
and the TM theories. We have 

5f(c) = l-2^M and S f(u,) = l-2gg4 (8) 



Df{u>) 



where Cf(u), Df{u), Cf(u) and Df(u) are 2 x 2 de- 
terminants which are explicitly given by 

Cf(oj) = ^ e(ui) / p,(uj)J' e (Lua/ c) Je (n(w) ua/c) 

— j£(uja/c)J' e (n(uj)uja/c), (9a) 

Df{oj) — e{uj) / p{uo)H^ (uja/c)Jg(n(uj)uja/c) 



—Hg {uaj c)Jg(n(u>) ujcl/c), 



(9b) 
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FIG. 1: TE theory, a) Total cross section cry . b) Scattering 
resonances in the complex uia/c plane. 
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FIG. 2: TM theory, a) Total cross section crfl. b) Scattering 
resonances in the complex um/c plane. 



and 



Cf(uj) = y fi(u))/e(u))J' e (u>a/c)Je(n(uj) uia/c) 

— J^(tva/c) J^(n(w) wa/c), (10a) 

Df(ui) = y/ fj,(u) I e(u))H^ (uja/c)j£(n(uj) ua/c) 

-HP(ua/c)J' e ( ti(lu) luo,/ c). (10b) 

For both polarizations, the unitarity of the S matrisi*, 
which expresses the energy conservation, and the reci- 
procity property^, which is associated with time-reversal 
invariance, can be easily verified by using elementary 
properties of Bessel functions. 

As far as the scattering amplitude f(cu, 9) and the total 
scattering cross section per unit length of the cylinder 



gt(lu) are concerned, they are, respectively, given by 



l_^ 7 ,[S^)-l]cos(£0). (ll) 



i=0 



where "fi is the Neumann factor (70 = 1 and for I 7^ 0, 
"/i = 2) and 



<tt(lu) 



/ 8tt 



■Im 



(e- < -/V(w,fl = 0)). (12) 



It should be noted that the definition (|12|) assumes the 
unitarity of the S matrix. It is obtained by using the 
optical theorem 1 —. From the expressions (|8|)- (llL)fl of the 
S matrices and by using expressions (fTT|> and (|T2*|l . we 

obtain the total cross sections (lu) and cr^(u>) for the 
TE and TM theories. We have 



a% (lu) = 



£=a 



and 



k(cu) 



E 



It 



Df(w) 



Df(«) 



(13) 



(14) 



Here, it should be noted that our expressions for the to- 
tal cross sections (lu) and cr£(ui) are different from 
the expressions given by Kuzmiak and Maradudinii. It 
seems to us that these authors have incorrectly applied 
the boundary conditions for the electric and magnetic 
fields at the cylinder surface. 

In Figs.^, and|2t, we display the total cross sections 
for the TE and the TM theories. They are both plotted as 
functions of the reduced frequency Lua/c. In the two fig- 
ures, rapid variations of sharp characteristic shapes can 
be observed. For the H polarization, this strongly fluctu- 
ating behavior is localized within and slightly around the 
frequency range where the cylinder presents left-handed 
behavior, while for the E polarization, it is totally lo- 
calized within that frequency range. For both polariza- 
tions, such strongly fluctuating behavior is due to the 
scattering resonances associated with the long-lived res- 
onant modes of the cylinder - i.e., the long-lived resonant 
states of the photon-cylinder system. These resonances 
are the poles of the S matrix lying in the fourth quadrant 
of the complex lu plane near the real lu axis. Resonances 
arc determined by solving 



Df(u)=0 for £eN 
for the TE theory and 

Df(u)=0 for feN 



(15) 



(16) 



for the TM theory. In Figs. ^ and|2b, resonances are 
exhibited for both theories. For certain frequencies, we 
can clearly observe a one-to-one correspondence between 
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the peaks of ctt(^) and the resonances near the real wa/c 
axis but in general the situation seems very confused. 
This is due to the profusion of long-lived resonant modes 
in and around the frequency range where the cylinder 
presents left-handed behavior. Furthermore, by zooming 
in on the distribution of resonances in regions close to the 
real axis of the complex to plane, we have also observed 
accumulations of resonances for large values of I: 

- For the TE theory, there exists an accumulation 
of resonances which converges to the limiting frequency 
lu s satisfying 



and given by 



e(u a ) + 1 = 



V2~ 



(17) 



(18) 



We have for the corresponding numerical reduced fre- 
quency uj s a/c w 7.806. 

- For the TM theory, there exists an accumulation 
of resonances at the limiting frequency u) / satisfying 



fj,(u f ) + 1 = 



and given by 



LOf = IVo 



2-F 



(19) 



(20) 



We have for the corresponding numerical reduced fre- 
quency Ufa/c 6.172. 

- For both theories, there exists an accumulation 
of resonances at the pole ojq of fj,(uj) and which corre- 
sponds more precisely to 



-co. 



(21) 



In summary as far as the spectrum of resonances is 
concerned, the left-handed cylinder is a physical sys- 
tem much richer than the metallic or the semiconducting 
cylinder (see Figs. 1 and 2 of Ref. 40 and the discussion at 
the end of Sec. II of that reference) and this is certainly 
very interesting for practical applications of left-handed 
electromagnetism. With this aim in view, it is necessary 
to understand the resonance spectrum of the left-handed 
cylinder from a physical point of view - i.e., to decode 
the underlying physics. That is what we shall do in the 
next sections. More precisely, we shall prove that this 
resonance spectrum is generated by SP's orbiting around 
the cylinder and we shall provide a numerical and a the- 
oretical description of these surface waves. 

It should be noted that the spectrum of resonances lies 
beyond the frequency range where n(uj) < (long-wave 
and short-wave scattering). The corresponding resonant 
modes arc associated with bulk polaritons. Their com- 
plex resonances appear in Figs, lb and 2b. Their imag- 
inary parts are much larger in modulus than those asso- 
ciated with the resonant modes generated by SP's. As 
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FIG. 3: The Watson integration contour. 



a consequence, they have a shorter lifetime, they do not 
play a significant role in the scattering process (see Fig. la 
and Fig. 2a beyond the frequency range where n(uj) < 0) 
and so they are much less interesting with in mind prac- 
tical applications. For all these reasons, we focus our 
interest only on the resonant modes associated with sur- 
face polaritons. 



III. REGGE POLES AND SURFACE 
POLARITONS 

Using the CAM method, we shall first provide a phys- 
ical picture of the scattering of electromagnetic waves by 
the left-handed cylinder in term of diffraction by surface 
waves. By means of a Watson transformation 47 applied 
to the scattering amplitude (|llfl . we can write 



2irk(Lu) 



V I £*M^ C0B[A(7T-fl)]rfA. 

sm7rA 



Here, in order to simplify the notation, we have not spec- 
ified the polarization. In Eq. I|22|l . C is the integration 
contour in the complex A plane (CAM plane) illustrated 
in Fig.|3|and which encircles the real axis in the clockwise 
sense. V which stands for Cauchy's principal value at the 
origin is used in order to reproduce the Neumann factor. 
The Watson transformation has permitted us to replace 
the ordinary angular momentum i by the complex angu- 
lar momentum A. S\(uj) is now an analytic extension of 
Se(u)) into the complex A plane which is regular in the 
vicinity of the positive real A axis. Using Cauchy's the- 
orem and by noting that inside the contour C the only 
singularities of the integrand in (|22|l are the integers, we 
can easily recover 111 Hi from Q22|). 

We can then deform the path of integration in i|22|) 
taking into account the possible singularities. The only 
singularities that are encountered are the poles of the 
S matrix lying in the CAM plane. They are known as 
Regge polesi2ii£ and are determined by solving 



Df (u) = for uj > 



(23) 
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FIG. 4: Regge poles in the CAM plane for the TE theory, 
a) The distribution corresponds to cua/c = 5.6. b) The dis- 
tribution corresponds to cua/c = 6.2. c) The distribution 
corresponds to cua/c = 7.0. 



FIG. 5: Regge poles in the CAM plane for the TM the- 
ory, a) The distribution corresponds to cua/c — 5.6. b) The 
distribution corresponds to cua/c — 6.2. c) The distribution 
corresponds to cua/c — 7.0. 



for the TE theory and 

Df{u) = for uj > (24) 

for the TM theory. Figures Q] and [5] exhibit the distribu- 
tion of Regge poles for both theories for three different 
reduced frequencies lying in the frequency region where 
ti(uj) < 0. At first sight, these Regge pole distributions 
are more complicated than the distributions associated 
with the metallic and the semiconducting cylinders stud- 
ied in Ref.^3- However, we have identified and indicated 
some particular Regge poles which, as we shall see be- 



low, are associated with surface waves orbiting around 
the cylinder and which explain its resonant behavior. For 
both polarizations: 

- One of these Regge poles is associated with the 
SP denoted SPoo which, as we shall show in Sec. V, 
corresponds, in the large-radius limit (i.e., for a — > oo), 
to a SP which is supported by the plane interface and 
which has been theoretically described in Refs. ISfilEvillssL 

- The other Regge poles are associated with an 
infinite family of SP's of whispering-gallery type denoted 
by WGSP n with n e N and which have no analogs in the 
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plane interface case (see Sec. V). 

For the TE theory, the Regge pole Asp^ lies in the first 
quadrant of the CAM plane while for the TM theory it 
lies in the fourth quadrant of that plane. For both theo- 
ries, the Regge poles Awgsp„ he in the fourth quadrant 
of the CAM plane but as u> — ► u>b they migrate to the 
third quadrant of that plane where they do not play any 
role. 

By Cauchy's theorem we can now extract from <|22[) 
the contribution of a residue series over Regge poles. For 
each Regge pole Asp(w) (here SP stands for SP^ as well 
for WGSP„), we capture an associated contribution given 
by 



/sp(w,0) 



' 2-7T rsp(w) 
ik(u>) sin [7rAsp(w)] 



cos [A S p(w)(tt 



where r SP {uj) = residue (SaM) a=Asp(w) . 
By using 



(25) 



+00 

sin it \ z — ' 



a +i7r(2m+l)A 



m—0 



which is true if Im A > 0, we can write 



/sp(w,0) = - 1 



' 2in 



r S p(w) 



+00 



x ^ e -HA S p(a>)(0+2m7r) + 



= +iA S p(w)(27r-6/+2m7r) 



(26) 



when the Regge pole Asp lies in the first quadrant of the 
CAM plane. Expression (|26|l is therefore associated with 
the surface wave SP00 of the TE theory. By using 



1 



sin 7rA 



+00 

= 2^ e- M(2m+1)A 

m=0 



which is true if Im A < 0, we can write 
f SP (u>, 6) = + ,1 --^P- r SP (u) 




x ^2 ( e ~ iX , ' ,ui ' ii( ' ,+2 '" :7 ' 

m=0 



-iAsp(w)(27r-e+2m7r) 



(27) 



when the Regge pole Asp lies in the fourth quadrant of 
the CAM plane. Expression l|27(l is therefore associated 
with the surface waves WGSP„ with n e N of the TE the- 
ory as well as with the surface waves SPoo and WGSP„ 
with n e N of the TM theory. 

In Eqs. I|26|) and l|27|) . exponential terms correspond 
to diffractive contributions. This clearly appears by tak- 
ing into account the time dependence exp(— iu>t). The 



physical interpretations slightly differ according to the 
position of the Regge pole A§p in the complex A plane: 

When the Regge pole Asp(w) lies in the 
first quadrant of the CAM plane - i.e., when we 
can use Eq. (|2"6*|) - the term exp[zAsp(w)(#)] (resp. 
exp[iAsp(w)(27r— 9)]) describes the SP propagating coun- 
terclockwise (resp. clockwise) around the cylinder and 
Re Asp (u>) represents its azimuthal propagation constant 
while Im Asp(^) is its damping constant. The cor- 
responding exponential decay reads exp[— Im Agp(w)#] 
(resp. exp[— Im Asp(w)(27r — 8)]). 

When the Regge pole Agp(w) lies in the 
fourth quadrant of the CAM plane - i.e., when we 
can use Eq. (|27J) - the term exp[— i\sp(ui)(6)} (resp. 
exp[— iAsp(w)(27r — 9)]) describes the SP propagating 
clockwise (resp. counterclockwise) around the cylinder 
and Re Asp(^) represents its azimuthal propagation con- 
stant while — Im Asp(w) is its damping constant. The 
corresponding exponential decay reads exp[+Im Asp(w)#] 
(resp. exp[+Im Asp (u>) (2n — 0)] ) . 

Finally, in both cases, the sum over m in l|2t)|) and l|27|) 
takes into account the multiple circumnavigations of the 
surface wave around the cylinder as well as the associated 
radiation damping. 

From the previous discussion, it is important to keep in 
mind that the function Re Asp(w) provides the dispersion 
relation for the SP associated with the Regge pole Asp 
and that the phase velocity v p and the group velocity v g 
of that SP are therefore given by 



ReAsp(w) 



and v g — 



d ReAsp(w) 



(28) 



Here we have taken into account the fact that the SP is 
supported by the cylinder surface at p = a and therefore 
that its wave number is given by 



fcsp(w) 



ReAsp(w) 



(29) 



Moreover, it is also important to note that the Regge 
poles Asp^ and Awgsp„ with n £ N are always close to 
the real axis in the complex A plane. As a consequence, 
they all correspond to SP's which are slightly attenu- 
ated during their propagation and which contribute sig- 
nificantly to the scattering process and to the resonance 
excitation mechanism. 

As ui varies, a given Regge pole Asp(w) describes a 
curve in the CAM plane. Such a curve is called a Regge 
trajectory^. From a physical point of view, the Regge 
trajectory corresponding to Asp(w) provides the disper- 
sion relation as well as the damping of the surface wave 
associated with this pole. In Figs. 16191 we have displayed 
the Regge trajectories of some SP's for the TE and TM 
theories. They have been obtained by solving numeri- 
cally Eqs. ill'- ill and l|24|) . We can observe some interesting 
features: 

- The dispersion curve for the surface wave SPoo 
of the TE theory is a positive and monotonically increas- 
ing function of lu. As a consequence, the associated group 
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FIG. 6: Regge trajectory for the Regge pole associated with 
SPoo (TE theory). 
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with the first three whispering- gallery SP's (TE theory). 



FIG. 9: Regge trajectories for the Regge poles associated 
with the first three whispering- gallery SP's (TM theory). 



and phase velocities given by Eq. (|28|) are both positive 
and SPoo has an ordinary behavior. It should be also 
noted that this SP exists in the frequency range ui €]0, lu s [ 
and therefore in the range ]wo 5 w b[ where the refraction 
index is negative but also outside this range. However, 
for low values of u>, its damping becomes very high and 
thus this surface wave has a negligible role in the scatter- 
ing process and in the resonance excitation mechanism. 
Furthermore, it should be noted that as u) — > u) s the dis- 
persion curve increases indefinitely. In the next section, 
this result will permit us to explain the accumulation of 



resonances which converges to the limiting frequency ui s 
for the TE theory. 

- The dispersion curve for the surface wave SPoo 
of the TM theory is a positive and monotonically de- 
creasing function of cj. As a consequence, the associ- 
ated phase velocity is positive while the group velocity 
is negative (see Eq. (J2HJ)- SP^ has a "left-handed be- 
havior" . It should be also noted that this SP only exists 
in the frequency range ]w/,o;fc[ which is included in the 
frequency range ujb[ where the refraction index is neg- 
ative, that its damping is always weak and thus that this 
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surface wave always plays a significant role in the scatter- 
ing process and in the resonance excitation mechanism. 
Finally, it should be noted that as w — > LOf the disper- 
sion curve increases indefinitely. In the next section, this 
result will permit us to explain the accumulation of reso- 
nances which converges to the limiting frequency uif for 
the TM theory. 

- As far as the surface waves WGSP„ with n £ N 
of the TE and TM theories are concerned, it seems at 
first sight they present a behavior which is rather inde- 
pendent of the polarization. It should be also noted that 
the real part of a given Regge pole Awgsp„ vanishes for 
a frequency in the frequency range ]cjo,w&[ and becomes 
negative. The Regge pole then lies in the third quad- 
rant of the CAM plane and is not taken into account 
by the theory previously developed. So we can consider 
that the surface waves WGSP„ with n £ N only exist 
in a sub-domain of the frequency range ]wo,^b[ where 
the refraction index is negative. The dispersion rela- 
tions of all these surface waves are positive but mono- 
tonically decreasing functions. Their group velocities are 
always negative while their phase velocities are positive 
(see Eq. pg)). All these SP's thus have a "left-handed 
behavior" . Furthermore, because the dampings of these 
surface waves are always weak, they all play a signifi- 
cant role in the scattering process and in the resonance 
excitation mechanism. Finally, it should be noted that 
as u — > wo, the dispersion curves increase indefinitely 
In the next section, this result will permit us to explain 
the accumulation of resonances which converges to the 
limiting frequency loq for the TE and TM theories. 



IV. FROM SURFACE POLARITONS TO 
RESONANCES 

From now on, we shall focus our attention on the physi- 
cal interpretation of the long-lived resonant modes whose 
excitation frequencies are localized within or slightly 
around the frequency range ujq < lu < LOb where the 
cylinder presents left-handed behavior. In this section, 
we shall prove that they are generated by the SP's de- 
scribed in Sec. Ill and for this reason we have called 
them RSPM's. 

We first consider that the solutions of l|15(l or l|16(l are 

denoted by loi p = — iTi p /2 where > and 

Tg p > 0, the index p permitting us to distinguish be- 
tween the different roots of (|15l) or l|l(jD for a given I. 
We also assume that in the immediate neighborhood of 
the resonance uji p , Sf(ui) or Sf(ui) and therefore the 
corresponding scattering amplitudes obtained from 
present a Breit-Wigner type behavior - i.e., are propor- 
tional to 



Tip/ 2 



(30) 



The resonant behavior of the cylinder (or more precisely 
of the cylinder-photon system) can then be understood 
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FIG. 10: Resonances generated by SPoo (TE theory). 
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in terms of SP Regge trajectories. Let us consider a 
given SP (we do not precise its polarization) described 
by the Regge pole denoted Asp(w). When the quantity 
Re Asp(w) coincides with an integer, a resonance occurs 
which is produced by a constructive interference between 
the different components of the surface wave, each com- 
ponent corresponding to a different number of circum- 
navigations of the cylinder (see Eqs. I|2vj|> and J23). Res- 
onance excitation frequencies w^gp associated with that 
SP are therefore obtained from the Bohr-Sommerfeld- 
type quantization condition 



ReA 



SP 



(J 0) \ 



= 0,1,2,. 



(31) 



Thus, they can be obtained numerically from the disper- 
sion curves. Now, by assuming that to is in the neighbor- 
hood of w^gp, we can expand Asp(w) in a Taylor series 



about ujpgp and write 



AspM ^1 + 



dReAsp(w) 



du) 



+«ImAsp(^g P ) + i 



dlmAsp(w) 



dio 



{co - uj esp ) 
(32) 



Then, by replacing (|32|) in the term sin [7tAsp(w)] of J2HJ, 
we can see that /sp(w, 69) presents a resonant behavior 
given by the Breit-Wigner formula (|30fl with 



P 



/SP 



ImAsp (w) (dReAgp (w) / dto) 



(dRe\ SP (Lj)/dLj) 2 + (dImA SP (w)/dw) 2 



(33) 

Our reasoning also predicts a shift in the excitation fre- 
quency w^gp solution of Eq. Ij31|l which is given by 

ImAsp {lo) (dlmAsp (w) / dto) 



A^sp - 



(dReX SP (uj)/duj) 2 + (dlm\ SP (Lo) /dio) 



In fact, such a shift is always negligible and it can be 
forgotten. Furthermore, it should be noted that in some 
frequency range we have 



\dRe\sp(uj)/duj\ ^> \dlm\sp (uj)/dt 
so that then reduces to 



r^sp 



Im Asp(w) 



d ReAsp(w)/dw 



,(0) 



(35) 



(36) 



Equations H31|) and (|33|l (or l|3()l) ) are two semiclassi- 
cal formulas which permit us to determine the location 
of the resonances from the Regge trajectory of Agp. In 
Figs. llOlTSI we present samples of complex frequencies for 
the RSPM's associated with the surface waves SP^ and 
WGSP„ with n = 1, 2 and 3. They have been calculated 
from the semiclassical formulas l|31fl and l|33|l by using 
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the Regge trajectories determined numerically by solv- 
ing Eqs. l(23l and (see Figs. I6I9|) . A comparison be- 
tween the semiclassical spectra and the "exact ones" (cal- 
culated by solving numerically Eqs. (fPo)) or ffT^i ) shows 
a very good agreement. Moreover, we can also observe 
some interesting features: 

- The resonance spectrum associated with the sur- 
face wave SPoo of the TE theory (see Fig. ITU|l extends 
beyond the frequency range ]w ,Wfc[ where the cylinder 
presents left-handed behavior because SPoo exists for 
u> G]0,w s [. Furthermore, inserted into the semiclassical 
formulas l|31l) and . the behavior of the Regge tra- 
jectory of Asp^ near lo s easily explains the existence of 
the family of resonances close to the real axis of the com- 
plex lo plane which converges for large I to the limiting 
frequency lo s . 

- The resonance spectrum associated with the sur- 
face wave SPoo of the TM theory (see Fig. El fully lies 
inside the frequency range Jc^Oj^bl where the cylinder 
presents left-handed behavior because SPoo exists only 
in that range. Furthermore, inserted into the semiclas- 
sical formulas 1)31(1 and J33J), the behavior of the Regge 
trajectory of Asp^ near LOf explains the existence of the 
family of resonances close to the real axis of the com- 
plex lo plane which converges for large t to the limiting 
frequency u> / . 

- The resonance spectra associated with the sur- 
face waves WGSP„ with neNof the TE and TM the- 
ories (see Figs. and II 3[) fully lie inside the frequency 
range ]ojq, LOb[ where the cylinder presents left-handed be- 
havior because all these surface waves exist only in that 
range. Furthermore, inserted into the semiclassical for- 
mulas (|31|l and (|33[) . the behavior of the Regge trajectory 
of a given Regge pole Awgsp„ near loq explains the ex- 
istence of a corresponding family of resonances close to 
the real axis of the complex lo plane which converges for 
large £ to the limiting frequency loq. Since there is an 
infinity of surface waves WGSP„, the accumulation of 
resonances in Fig. ^ (resp. Fig. Eb) is more intense at 
the point loq than at the point lo s (resp. ujf). 

In conclusion, we have established a connection be- 
tween the complex frequencies of the long-lived resonant 
modes (or RSPM's) of the cylinder and the SP's noted 
SPoo and WGSP„ with n 6 N which are supported by 
its surface. In other words, in spite of the great confu- 
sion which seems to prevail in the resonance spectrum 
of the left-handed cylinder (see Sec. II), we have been 
able to fully classify and physically interpret the reso- 
nances thanks to CAM techniques. We now invite the 
reader to look at Figs. 1 141 and 1 151 where we have zoomed 
in on Figs.^i and|2Ji. On the total cross sections er^ (lo) 
and cr|? (lo) we have identified the peaks corresponding to 
resonances and, for each one, we have specified the SP 
which has generated it as well as the associated "quan- 
tum number" £. 



Cylinder TE 




FIG. 14: Zoom-in on the total cross section <x<f? . 




FIG. 15: Zoom-in on the total cross section a T . 



V. ASYMPTOTICS FOR SURFACE 
POLARITONS AND PHYSICAL DESCRIPTION 

In order to obtain a deeper physical understanding of 
the SP's orbiting around the left-handed cylinder and to 
justify the terminology previously used, we must "ana- 
lytically" solve Eqs. (|33J) and l|23} for A = Asp or equiv- 
alently 



J' XBP (\n(uj)\coa/c) _ U(lo) H^ooa/c) 



J Xsp (\n(Lo)\Loa/c) 
for the TE theory and 

^A SP (K W )I wa / c ) 



MM HQ (uoa/c) 



l^uo)H^l{Loa/c) 



J\sA\n(u)\ua/c) y e(uo) ( wa / c ) 



(37) 



(38) 



for the TM theory Here, Eqs. J57JI and iJSSJl have been 
obtained from Eqs. (|23() and l|24ll by using the rela- 
tions Jx(-z) = e i7rX J x (z) and J' x (-z) = -e" x J' x (z) (sec 
Ref. l46|) . These two equations can be solved approxi- 
matively by using asymptotic analysis. More precisely, 
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by considering adequate asymptotic expansions for the 
Bessel functions occurring in (|37[) and (J3SJ, we shall be 
able to provide analytic expressions for the Regge poles 
Asp^ {lo) and A\vgsp„ (<^) for the two polarizations. Of 
course, the main difficulty with this method concerns the 
choice of the asymptotic expansions for the Bessel func- 
tions. Such a choice strongly depends on the relative 
positions of the arguments \n{oj)\u>a/ c and toa/c with re- 
spect to the complex order Asp- In order to simplify the 
discussion, we choose to describe theoretically the SP's 
in the frequency ranges where they generate the RSPM's 
with the longest lifetime (such modes are the most im- 
portant from the physical point of view). In other words, 
we shall seek Asp^ (w) of the TE theory for oj in the 
neighborhood of ui s , Asp^w) of the TM theory for u> in 
the neighborhood of Uf and A\yGSP„(kO of the TE and 
TM theories for ui in the neighborhood of ujq ■ In fact, 
in spite of these restrictions, we shall obtain asymptotic 
results valid in large frequency ranges. 

Let us first consider the Regge pole associated with 
SPoo of the TE theory. We assume ui in the neighborhood 
of uj s and then we can also assume that ReAsp oo (w) > 
Lua/c > \n(aj)\uja/c and formally that ReAsp oo (o;) ^> 1 
and ReAsp =c (w) ^> ImAsp 00 (w). The configuration there 
is displayed in Fig.^J:. As a consequence, we can use the 
Debye asymptotic expansions for J\(z) and H^\z) valid 
for large orders (see Appendix A of Ref. or Ref. E^) 
in the form 

J x {z)~ l -A{\,z)e +a{x > z \ (39) 
f\z)~-tA(\,z)e- a ^\ (40) 



n Im A 



H 



where 

A(X,z) 



1/2 



(A 2 - z 2 r 1/4 , (41a) 

v2 _ 



k(A, z) = (A 2 - z 2 ) 1 ^ _ A In (^ A+(A2 ~ Z2)1 2 ) 



(41b) 



Thus, we have 

J A SPoo (l n Hk a /c) 
■hs?^(\n{uj)\uja/c) 



[Al Poo -(|nH| W a/ C ) 2 ] 1/2 



and 



H 



(1)' 
a sp 



(ua/c) 



[A 



2 

SP. = 



\n{<jj)\(uia/ c) 



{(jja/cY 



<L>«/ C ) 



(u>a/c) 



(42) 



(43) 



Now, by inserting (g^J and l|4H|l into Eq. IpTfjl , we obtain 
an equation which can be easily solved and we deduce 



AsPooH 



e(u)\(\e(L,)\+»(u)) 
e 2 (uj) - 1 



(44) 
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FIG. 16: The relative positions, in the A complex plane, of 
the reduced frequency f3 — uja/c and the Regge pole Asp^ of 
the TE theory. 



We have obtained an asymptotic expansion for 
AsPoo or, more exactly, for the real part of that func- 
tion. Indeed, it should be noted that the right-hand 
side of is purely real. The perturbative method 
previously used did not permit us to extract the small 
imaginary part of Asp^ (to). Even by taking into account 
higher orders in the asymptotic expansions (|39|) and 14UII . 
we are not able to provide the imaginary part of Asp^ (w) 
because it corresponds to an exponentially small contri- 
bution which lies beyond all orders in perturbation theory 
and which can be captured only by carefully taking into 
account the Stokes phenomenon^2^i. We have already 
encountered such a problem in Ref. and, mutatis mu- 
tandis, we shall solve it in the same way. Instead of l|4(J [) . 

we must use the Debye asymptotic expansion of (z) 
in the form 



h£\z)~ - iA(\, z)e- a{x ' z) (l + ...) 



+S[a{X,z)]A(X,z)(l 



)e 



a(\,z) 



(45) 



In the right-hand side of 14511 . the first term is the usual 
Debye asymptotic expansion truncated near its least 
term. The second one is obtained by decoding the diver- 
gent tail o f that asymptotic expansion. This can be done 
(see Refs. I42l43|) by Borel summation after exploiting a 
resurgence formula discovered by Dingle^i. In the region 
of the A complex plane where the Regge pole Asp^ lies 
(see Fig. EJ|, we have Re a < 0. As a consequence, the 
first term of the right-hand side of l|45|l is the dominant 
contribution while the second one is a subdominant term 
which can be forgotten when |A| — > oo. That is what we 
did previously by using (|40|l . The Stokes multiplier func- 
tion S [a(X, z)] is a complicated function involving the 
exponential integral function E\. It goes continuously 
from to 1 at the crossing of the Stokes line Im a — 
emerging from the turning point z = (3 (see Fig. 116(1 . Be- 
low the Stokes line, it rapidly vanishes. On the Stokes 
line it is equal to 1/2 and above the Stokes line it rapidly 
becomes equal to 1. It thus describes the rapid but con- 



13 



tinuous birth of the subdominant contribution near the 
Stokes line^Si. From (|45|) we can now write 



30 r 



H 



(i)' 

Asp 



(ioa/c) [A 



sp 



(ua/c) 



Hi 1 ! {uoa/c) 



^SF 



X 1 



(coa/c) 

2iS [a(X SPoo , toa/c)] e 2a( - Xs ^ 

(46) 

instead of ((43(1 . Because we can assume that Asp^fy) is 
very close to the Stokes line emerging from the turning 
point uja/c, we can consider that S [a(\gp x (u/), ua/cj] = 
1/2. Then, by inserting and (gHJl into Eq. l(37|l. we 
obtain 




ReA 



ImAsp^ (uj) 



eH|(|e( W )|+/iM) 



e 2 (w) - 1 



(47a) 



e 2 H - 1 



^ (ReA SPoo M) -(wa/c) 

ReA SPoo (w) 
x exp[2a(ReAsp <X) (w), ua/c)]. 



(47b) 



Equations (|47a|l and (|47b|) provide analytic expressions 
for the dispersion relation and the damping of the surface 
polariton SPoo of the TE theory. The following impor- 
tant features must be noted: 

- The wave number ksp^ {u>) associated with this 
SP is obtained from l|47a|l and ((29(1 and is given by 



fcsp^ (w) 




e( W )|(|eH|+A*H) 



e 2 (w) - 1 



(48) 



This expressio n is the usual dispersion relation found in 
Refs. M.37 38 for the p-polarized SP - i.e., the SP for 
which the magnetic field H is normal to the incidence 
plane - supported by the flat interface. Here, we have 
recovered exactly the same dispersion relation because 
we have limited the perturbative resolution of Eq. I|37|) 
to the lowest order. By taking into account higher orders 
in the asymptotic expansions 1)39(1 and 145(1 . we could 
obtain corrections for ((481) which vanish for a — > oo - 
i.e., in the flat interface limit. The surface wave of the 
TE theory denoted by SP^ is therefore the counterpart 
of the p-polarized SP supported by the flat interface. 

- The imaginary part l(47b(l of Asp^ vanishes for 
a — > oo — i.e., in the flat interface limit. Then, SP^ has 
no damping like the p-polarized SP supported by the flat 
interface. This result strengthens our previous claim. 

- The function ReAsp OQ (w) given by l(47a|) has a 
simple pole when e(uj) + 1 = i.e., for lu = u> s . Fur- 
thermore, the imaginary part l(47b|l of Asp^ vanishes for 
w = lu s . These two results justify all our previous re- 
marks concerning the accumulation of resonances which 
converges to the limiting frequency lo s . 
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- We have numerically tested formulas ((47a(l and 
((47b(l (see Fig. I17|) . They provide rather good approxi- 
mations for ReAsp^ (u>) in a large frequency range and for 
ImAsp^ (tv) in the neighborhood of ui s . It would be possi- 
ble to improve ((47b|) by taking into account the variation 
of the Stokes multiplier and to obtain a formula valid 
in a larger frequency domain but then the expression of 
ImAsp^ (tv) would become much more complicated. 
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Let us now consider the Regge pole associated with 
SPoo of the TM theory. We must solve Eq. (gSJ which 
only differs from Eq. (|37|l by the factor in front of the 
Bessel functions. We assume u) in the neighborhood of 
ujf and then we can also assume that ReAsp oo (aj) > 
\n(uj)\uja/c > ua/c and formally that ReAsp DO (w) 3> 1 
and ReAsp =c (w) 3> ImAsp 00 (w). The configuration there 
is displayed in Fig. |Sp. As a consequence, we can use 
again the Debye asymptotic expansions for J\(z) and 

H^\z) given in i|39[l and (|45|l and the resolution of 
Eq. (2H1) can be modelled on that of Eq. l(37|). Formu- 
las (l4*2l and (|4l)|> are still valid and the Stokes multiplier 
in (I46() can be taken equal to 1/2. Indeed, in this case 
Asp^ (lo) is below the Stokes line emerging from the turn- 
ing point uja/c but it remains very close to this Stokes 
line. Then, by inserting g2l an d GUI mto Ec i- EH), we 
obtain 



A lm A 



ReA SI 



( w ) ~ (—) / l^MKI/^HI + e H) 
V c / y m 2 ( w ) — i 



(49a) 



ImA 



SP: 



M 2 H - 1 



ReAsp^ (w) 
x exp[2a(ReAsp OQ (w), wo/c)]. 



(49b) 



Equations (|49a|) and 149b|l provide analytic expressions 
for the dispersion relation and the damping of the surface 
polariton SPoo of the TM theory. The following impor- 
tant features must be also noted: 

- The wave number fcsp^ (lo) associated with this 
SP is obtained from 149a(l and (|29|l and is given by 



^SPooM 




mMKImMI + e(w)) 



m 2 H - 1 



(50) 



This expressio n is the usual dispersion relation found in 
Refs. I3fil37l38l for the s-polarized SP - i.e., the SP for 
which the electric field E is normal to the incidence plane 
- supported by the flat interface. The surface wave of the 
TM theory denoted by SPoo is therefore the counterpart 
of the s-polarized SP supported by the flat interface. 

- The imaginary part l)49bj) of Asp^ vanishes for 
a — > oo - i.e. in the flat interface limit. Then, SPoo has 
no damping like the s-polarized SP supported by the flat 
interface. This result strengthens our previous claim. 

- The function ReAsp^ (u>) given by l)49al) has a 
simple pole when + 1 = 0- i.e., for lu = ujf. Fur- 
thermore, the imaginary part l|49b(l of Asp^ vanishes for 
lo = u>f. These two results justify all our previous re- 
marks concerning the accumulation of resonances which 
converges to the limiting frequency tOf. 

- We have numerically tested formulas (|49a|l and 
(|49b(l (see Fig. I18fl . They provide good approximations 
for ReAsp OQ (w) in the frequency range and for 
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FIG. 19: The relative positions, in the A complex plane, of 
the reduced frequencies (5 — uia/c and 7 = \n(u))\u>a/ 'c and of 
the Regge pole Awgspi of the TE and TM theories. 



ImAsp 00 (w) in the neighborhood of tOf. Of course, it 
would be possible to improve l|49bl) by taking into account 
the variation of the Stokes multiplier and to obtain a 
formula for the imaginary part of Agp^ (to) valid in a 
larger frequency domain. 

Let us finally consider the Regge poles associated 
with the surface waves WGSP„ for the TE and TM 
theories. We must now solve Eq. Ij37|l and Eq. i|38|) 
for Asp = Awgsp„ by assuming lo in the neighbor- 
hood of <x>o- We are in the configurations described in 
Fig. Hi for the TE theory and Fig. 0} for the TM the- 
ory (see also Fig. I19f) . As far as the asymptotic ex- 
pansion for H\(z) is concerned, we can still use l|40[) 
as well as its modification (I45() because we can as- 
sume that ReAwGSP„ (w) > wa/c and formally that 
ReA W GSP„M > 1 and ReA W GSP„M > ImA W GSP„M- 
We then obtain 



H 



(i)' 

AwGSPt 



(cua/c) 



[A 2 



WGSP„ 



{uia/cY 



H 



1wosp>«/ C ) ( " a/c) 

x (l - 2iS [a(A W GSP„ , wa/ c)] . 



,2a(A W GSP„ ,i^a/c) 



(51) 



As far as the asymptotic expansion for J\{z) is concerned, 
the situation is a little bit more complicated. Awgsp„ is 
in the immediate neighborhood of \n(uj)\uja/c or more 
precisely in the Airy circle centered on \n(uj)\u>a/c (see 
Fig. I19fl . This is true for the Regge pole corresponding 
to n = 1 but it is also formally true for all the other 
Regge poles with n 7^ 1 when lo — > loq. Then, the Debye 
asymptotic expansion l|40l) for J\(z) fails because A ps z 
and we must use the uniform asymptotic expansion (see 
Appendix A of Ref. |H or Ref . 



J x (z)~(2/z) 1/3 Ai (2/z) 1/J (A-z) 



(52) 
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where Ai(z) denotes the Airy function^. Thus we have 



J 



AwGSP, 



(\n(uj)\u>a/c) 



J\wGsr n {\ n ( u )\ua/c) 



(2/(\n(u,)\u,a/c)) 1/3 



Ai 



1 /3 

(2/(\n(w)\wa/c)) (A W gsp„ - \n(u))pa/c) 



1 /3 

Ai (2/(|n(w)|wa/c)) (A W GSP n - \n(uj)\u)a/c) 

(53) 

We can now insert and J5TJ into Eqs. (|37|) and 
In order to solve the resulting equations, we first note 
that their rhs tend to infinity when lo — > luq. Hence, 
the argument of the Airy function is very close to the 
zeros x n of that function (let us recall here that the first 
three ones are xi w —2.3381..., x 2 ~ —4.0879... and x% « 
—5.5205...) and we then have 



Awgsp„M ~ \n(uj)\uja/c + 



\n(u)\u>a j c\ 1 



(54) 



Of course, Eq. (|54|) provides an approximation which can 
be easily improved by taking into account the subdom- 
inant contribution in (|51|) . We assume that the Stokes 
multiplier is equal to 1/2. This is true because the Regge 
pole Awgsp„ (w) is below but very close to the Stokes line 
emerging from the turning point Loa/c. We then obtain 



ReAwGSP„ (w) ~ \n(id)\u)a/c + 

for the TE and TM theories, 
ImAwGSP„ M >■ 



\n(u>)\(jja/ c\ 1/ ' 3 



(55a) 



e(tu) ( \n(io)\ oja/c\ 2 ^ 3 



(Re A W gsp„M) 2 - (uja/cf 



1/2 



(u;a/c) x n 

x exp[2a(Re A WG sp„ M, wa/c)] (55b) 
for the TE theory and 



ImA W GSP„ M 



/i(w) ^|n(cj)| uja/c^ 2 ^ 3 



(Re A WG sp„M) 2 ~ (uja/cY 



1/2 



(u;a/c) x n 

x exp [2a(Re A WG sp„ M, ua/c)] (55c) 

for the TM theory . 

Equations l|55a|) - (|55c|) provide analytic expressions for 
the dispersion relation and the damping of the surface 
polaritons WGSP n of the TE and TM theories. The cor- 
responding wave numbers /cwgsp„ (w) which are obtained 
from (|55a|l and (|29|) are given by 



fcwGSP„M ~ \n(uj)\(j/c + 



\ti{uj)\uj / c 
2^2 



1/3 



(56) 
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It should be noted that (|54|l and therefore l|55a|l have 
been obtained by following, mutatis mutandis, the cal- 
culations done by Rayleigh a long time ago in order 
to describe mathematically the whisperingjrallery phe- 



nomenon in acoustic 



,52,53 



(see also Ref. .54). We en- 
counter the same phenomenon but in the context of left- 
handed electromagnetism and the SP's described by the 
Regge poles Awgsp„ are surface waves of whispering- 
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gallery type. Here it should be recalled that a whispering- 
gallery surface wave propagates only close to a curved 
interface and mainly in the internal medium in a layer 
whose thickness is proportional to fc~ 2 / 3 with k its wave 
number. From a physical point of view, the surface po- 
laritons WGSP„ are therefore very different from the sur- 
face polaritons SPoo which exist on both sides of the in- 
terface. Moreover, the following important features must 
be also noted: 

- These SP's have no counterparts in the plane 
interface case. 

- The function ReA\yGSP„ (w) given by (|49a|l has 
a pole which is that of n{u) and therefore which corre- 
sponds to w = loq. Furthermore, the imaginary parts 
(I55bll and l|55c|) of Awgsp„ vanish for lu — oj a . These 
results justify all our previous remarks concerning the 
accumulations of resonances which converge to the lim- 
iting frequency loq. 

- We have numerically tested formulas l|55a|) - (|55c|) 
(see Figs . l2Dlandl?Tlfor WGSPi ) . They provide very good 
approximations for ReAwGSPi (w) in the full frequency 
range LUb[ where the cylinder presents left-handed be- 
havior. They also provide very good approximations for 
ImAwGSPi (w) in a rather large frequency range above the 
limiting frequency ojq. 



VI. CONCLUSION AND PERSPECTIVES 

In the present article, we have considered the scatter- 
ing of TE and TM electromagnetic waves by a cylinder 
fabricated from a left-handed material. We have mainly 
emphasized the resonant aspects of the problem. We 
have shown that the long-lived resonant modes can be 
classified into distinct families, each family being gener- 
ated by one SP and we have physically described all the 
SP's orbiting around the cylinder by providing, for each 
one, a numerical and a semiclassical description of its 
dispersion relation and its damping. 

We have also shown that the left-handed cylindrical 
interface can support both TE- and TM-polarized SP's. 
For each polarization, we have found that there exists 
a particular SP which corresponds, in the large-radius 
limit, to the SP which is supported by the plane in- 
terface and which has been theoretically described in 
Refs. 13603711381 But there also exists, for each polariza- 
tion, an infinite family of SP's of whispering-gallery type 
and these have no analogs in the plane interface case. 
The existence of these SP's of whispering-gallery type is 
certainly the main result of our article. Here, in order 
to appreciate it, it should be recalled that on a cylindri- 
cal metal-dielectric or semiconductor-dielectric interface, 
the SP's of whispering-gallery type are absent. There 
exists only one SP— . It can be excited in the H polar- 
ization configuration - i.e., by TE waves - and in the 
large-radius limit it corresponds to the SP which is sup- 
ported by the plane interface. From the point of view 
of SP physics, left-handed cylinders are therefore much 



richer systems than metallic or semiconducting cylinders 
and this could have important consequences in term of 
practical applications. 

We are now completing a theoretical and numerical 
study in order to provide a better physical understand- 
ing of the resonant modes of left-handed cylinders^. We 
more particularly focus on the repartition of the electro- 
magnetic energy density - both inside and outside the 
cylinder - associated with all the resonant modes. That 
enables us to physically distinguish on the one hand the 
resonant modes generated by bulk polaritons from those 
generated by SP's, and on the other hand, the resonant 
modes generated by SP's of SPoo type from those gener- 
ated by SP's of whispering-gallery type. 

It should be also noted that in our paper we have as- 
sumed for the effective electric permittivity and the ef- 
fective magnetic permeability of the left-handed material 
the simple expressions respectively given by © and 10} . 
In fact, actual left-handed materials have a more com- 
plicated permittivity. Indeed, they consist of periodic 
arrangements of split ring resonators and wires and it 
has been recently shown that the split ring resonators 
present a strong electric response in addition to the mag- 
netic one (see Ref. 56). In fact, it is possible to take 
into account this response by assuming that the effective 
electric permittivity is given by 56 



with wo 7^ while the effective magnetic permeability is 
still given by Such a modification does not change 
radically our previous analysis. Mutatis mutandis, our 
results remain valid. Indeed, we have theoretically - i.e., 
semiclassically - and numerically reexamined, in this new 
context, the resonant aspects of a left-handed cylinder 
linked with SP's. We have obtained the following results: 

- The surface wave SP,^ of the TE theory still ex- 
ists. It generates a family of resonances which converges 
to the limiting frequency oj s satisfying e(uj s ) + 1 = 
and which is now given by lo s — (lo v 2 + u>q)/2 instead 
ofgHJ- This SP is still described by Eqs. (|4Ta)l and 
(|47b(l with now e(ui) given by 1)57(1. Of course, its be- 
havior depends on the relative order of the characteristic 
frequencies ujq, loq, u>i,, oj p . In general, it presents a left- 
handed behavior if lu s lies in the frequency range where 
n(uj) < 0. 

- The surface wave SP^ of the TM theory still ex- 
ists. It generates a family of resonances which converges 
to the limiting frequency Uf satisfying /J-(ujf) + 1 = and 
which is still given by i|2(J|) . This SP is still described by 
Eqs. (|49a|l and (|49b(l and, in general, it presents a left- 
handed behavior if ujf lies in the frequency range where 
n(w) < 0. 

- For both polarizations, there still exists the 
family of SP's of whispering-gallery type described by 
Eqs. I(55a(l - ((55c(l in the neighborhood of the pole ujq of 
/i(w). They generate an infinity of resonances which con- 
verges to the limiting frequency luq. In general, they 
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present a left-handed behavior if loq lies in the frequency 
range where n(w) < 0. 

- But there also exists, for both polarizations, 
a new family of SP's of whispering-gallery type in the 
neighborhood of the pole Cjq of e(ui). They are also de- 
scribed by (|55a|l - (|55cjl (in this frequency range) and they 
generate an infinity of resonances which converges to the 
limiting frequency Cjq. In general, they present a left- 
handed behavior if ujq lies in the frequency range where 
n(cj) < 0. 

Finally, it is worthwhile pointing out the possible ap- 
plication of the CAM method in the context of multi- 
ple scattering by left-handed cylinders. This could be 



achieved by extending to electromagnetism of dispersive 
media the approaches previously developed in the context 
of quantum mechanics and acoustics (see, for example, 
Ref. [57 58). SP's probably induce a strong coupling be- 
tween the different cylinders and the CAM method could 
permit us to understand it and to quantify its intensity. 
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